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Abstract. We present B,V,I CCD photometry of about 
40000 stars in four regions of the Fornax dwarf spheroidal 
galaxy down to V ~ 23.5, the largest three-color data set ob- 
tained for this galaxy until now. The resultant color-magnitude 
diagrams, based on a wide color baseline, show a variety of fea- 
tures tracing the history of star formation of this dwarf galaxy. 
One of the most distinctive features in our diagrams is the 
conspicuous young main sequence, indicating recent star for- 
mation until approximately 2 x 10 8 yr ago. A plume of stars 
brighter than the red HB clump, with (B — I) ~ 0.5, trace the 
helium-burning phase of the young population. A comparison 
of the color and extension of this feature with model isochrones 
suggests a relatively metal-rich population ([Fe/H]~ —0.7) 
with age 300-400 Myr. This represents an important constraint 
for understanding the chemical enrichment history of Fornax. 
An extended upper AGB tail and a prominent red HB clump 
sign the presence of the well-known dominant intermediate- 
age population with an age range 2-10 Gyr, for which we have 
estimated a mean age 5.4 ± 1.7. About 0.2 mag below the red 
clump, an extended HB is indicative of an old population. We 
show that blue HB stars may be present in the outer regions. 
Together with previous detection of RR Lyrae, this provides 
evidence for a minority field population that is as old and metal- 
poor as that in the Fornax globular clusters. We have identified 
the AGB bump, a clustering of stars that occurs at the beginning 
of helium shell-burning evolution, at a luminosity My ~ —0.4. 
This is an example of the short-lived evolutionary phases that 
can be revealed in stellar populations using adequately large 
star data samples, whose measurements provide powerful tests 
of theoretical models. 

Based on precise detection of the tip of the RGB in a se- 
lected RGB sample, we measure a corrected distance modulus 
(m - M)o = 20.70 ± 0.12. An independent estimate of the 
distance to Fornax was also obtained from the mean magni- 
tude of old horizontal branch stars, yielding a distance modulus 
(to— M)q = 20.76±0.04, in good agreement with the distance 
estimated from the red giant branch tip and previous results. 
The large baseline of the (B — I) colors together with the size 



of the stellar sample allowed us to analyze in detail the color 
distribution of the red giant stars. We find that it can be approx- 
imately described as the superposition of two populations. The 
dominant component, comprising ~ 70% of the red giant stars, 
consists of relatively metal-enriched intermediate-age stars. Its 
mean metallicity is [Fe/H]=— 1.39 ± 0.15, based on a compar- 
ison of the fiducial locus of the bulk of the Fornax red giants 
with the homog eneou s Galactic globular cluster set of Da Costa 
&Armandroff( |l990| ). Once the younger mean age of Fornax is 
taken into account, our best estimate for the mean abundance of 
the bulk of the galaxy is [Fe/H]« —1.0 ± 0.15. The dominant 
intermediate-age component has an intrinsic color dispersion 
<Jo(B — I) = 0.06 ± 0.01 mag, corresponding to a relatively 
low abundance dispersion, crpe/H] = 0.12 ± 0.02 dex. Fur- 
ther, there is a distinct small population of red giants on the 
blue side of the RGB. While these stars could be either old or 
young red giants, we show that their spatial distribution is con- 
sistent with the radial gradient of old horizontal branch stars, 
and completely different from that of the younger population. 
This unambiguously qualifies them as old and metal-poor. This 
result clarifies the nature of the red giant branch of Fornax, 
suggesting that its exceptional color width is due to the pres- 
ence of two main populations yielding a large abundance range 
(—2.0 < [Fe/H] < —0.7). This evidence suggests a scenario 
in which the Fornax dSph started forming a stellar halo and its 
surrounding clusters together about 10-13 Gyr ago, followed 
by a major star formation epoch (probably with a discontinu- 
ous rate) after several Gyr. 

Key words: Galaxies: fundamental parameters - Galaxies: in- 
dividual: Fornax - Local Group - Galaxies: stellar content - 
Galaxies: structure 



1. Introduction 

An increasingly large number of investigations has recognized 
the importance of dwarf spheroidal galaxies for our under- 
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standing of galaxy formation and evolution (see Mateo 1998 
and Da Costa 1998| for recent reviews). While new studies of 



the central, densest regions of the more distant Local Group 
galaxies have benefited from the exceptional resolution of HST, 
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the nearby dwarf spheroidal satellites of the Milky Way can still 
be investigated in great detail using ground based wide-field 
data. Thanks to recent improvements in detector efficiency (in 
particular, in the blue part of the optical spectrum) and size, 
CMD's with long color baselines and high statistical signifi- 
cance can be obtained in affordable exposure times. 

The Fornax dwarf spheroidal ( dSph) galaxy represents one 
of the most interesting cases for studying the complexity of 
stellar populations in dwarf galaxies. This galaxy was one of 
the first dSph in which an intermediate age population was de- 
tected. The presence of upper asymptotic giant branch (AGB) 
stars, brighter and redder than the tip of the red giant branch 
(RGB), indicated that about 20% of the galaxy could be of in- 
termediate age (2 to 8 Gyr) (Aaronson & Mould [l980[ |l 98 5| ). 
Surveys for AGB stars led to discovery of 1 1 1 carbon stars, 
for which follow up near-infrared photometry indicated a wide 
range of bolometric luminosities, a mass dispersion among the 



turnoff was located at My ~ —1.4, implying a minimum age 
of ~ 10 8 yr. This young population is best shown by the recent 
wide-area survey of SHS98. Notwithstanding this young stel- 
lar component, Fornax appears to be devoid of any interstellar 
medium (ISM). A large-area search for neutral hydrogen has 



given no detectable Hi emission or absorption (Young 1999), 
the upper limit for Hi emission being 5 x 10 18 cm" 2 at the 
galaxy center. Thus the interstellar medium that must have been 
present a few 10 s yr ago to form stars, appears to be all gone. 
There is also the possibility that the ISM has been ionized and 
heated up by the interstellar UV field. However, this hypothe- 
sis conflicts with the lack of det ection of X-ray emission in the 
direction of Fornax (Gizis et al. 



1993) 



progenitors, and hence an age spread (Frogel et al. 1982; West- 
erlund et al. |1987| ; Lundgren [l990| ; Azzopardi et al. |1999| ). For- 
nax is also known to contain a planetary nebula whose abun- 
dance patterns are consistent with an origin from a second or 
third generation star (Danziger et al. 1978; Maran et al. 1984). 



The various stellar subpopulations in Fornax have different 
spatial distributions, which have been carefully investigated by 
SHS98. The oldest population, represented by the RR Lyrae 
variables, has the most extended distribution. The intermediate- 
age stars (red clump stars) are more centrally concentrated, 
while the young population of blue MS stars, as well the red- 
dest AGB (carbon) stars, are even more concentrated in a bar- 
like distribution roughly aligned in the EW direction, with the 
brightest stars located at the ends of the bar. Also the red clump 



The presence of such an intermediate-age population is con- population displays an asymmet rical structure (cf. Hodge |l96T 



Stetson et al. 1998). Most recently, an HST study of a cen 



firmed by a conspicuous red HB clump (Demers et al. |199j Eskridge p§| Demers et al rf^a wkh a peculiar « crescent - 



tral Fornax field sampling the main-sequence turnoffs of the 
intermediate-age and o ld populations has been carried out by 
Buonanno et al. ( 1999 ). The analysis of the resulting CMD has 
shown evidence for a star formation starting about 12 Gyr ago 
and continuing until 0.5 Gyr ago. A variable star formation rate 
is revealed by gaps between separate subgiant branches, and 
major star formation episodes probably occurred nearly 2.5, 4, 
and 7 Gyr ago. 

Also, Fornax certainly harbors an old stellar population, 
since it contains five globular clusters whose ages do not differ 



from those of M68 and M92 (Buonanno et al. 1998; Smith et al 



1998 ), except perhaps for cluster 4 that appears to be 2-3 Gyr 
younger (Buonanno et al. 1999J ). These clusters have unusually 



red horizontal branches for their low metallicity, with no coun- 
terparts in the outer Galactic halo or the Magellanic Clouds. 
Also for cluster 4, the r ecent WFPC2 color-magnitude dia- 
grams of Buonanno et al. ( 1999 ) unambiguously indicate a low 
metallicity, [Fe/H]« —2, although integrated spectra pointed 
to a m etallicity similar to that of field stars (Beauchamp et al. 
1995). An old population is present among the field stars of 



Fornax as well, as demonstrated by detection of a red horizon- 
tal branch slightly fainter t han th e red clump, a nd of RR Lyrae 
variables (Buonanno et al. 1985; Stetson et al. 1998, hereafter 



SHS98). 

Fornax also hosts a significant population of young stars. 



Buonanno et al. (1985) had already noticed a handful of faint 



blue stars in their CMD of the Fornax field, tentatively explained 
as belonging to a ~ 2 x 10 9 yr population (cf. also Gratton et 



al. 1986). The deeper CMD of Beauchamp et al. (1995) clearly 



revealed a young main-sequence, and comparison with theoret- 
ical isochrones indicated recent star formation. The brightest 



shape (SHS98). 

Despite all these pieces of knowledge accumulated in re- 
cent years, the star formation history of Fornax is not yet fully 
understood. Several questions need to be answered before a 
reconstruction of the star formation and chemical enrichment 
history of Fornax can be attempted. The metallicity should be 
measured for stellar populations of different age and location 
within the galaxy, so as to determine the run of metal enrich- 
ment as a function of time. The star formation history needs 
to be evaluated using critical features in the CMD as tracers of 
star formation at different epochs, to understand to what extent 
star formation proceeded continuously or in bursts, and how it 
propagated throughout the body of the galaxy. The nature of the 
wide red giant branch (RGB) is still quite puzzling, though all 
previous investigations agree on the fact that it is broader than 
expected on the basis of the photometric errors. Further, there 
is a lack of observational data with which to study features such 
as the RGB and AGB bumps or the precise location of central 
helium-burning stars as a function of age and metal abundance, 
as a test of stellar evolution models. Large field observations of 
Local Group dSph galaxies, being able to sample a significant 
number of stars, can address these issues. 

With these open questions in mind, we have investigated 
the stellar populations of Fornax as part of a wide-field study 
of nearby dwarf spheroidals. We present here a large area BVI 
photometric study of the Fornax field, yielding magnitudes and 
colors for about 40,000 stars down to ~ 2 mag below the hor- 
izontal branch, in four regions located at different distances 
from the galaxy center. The use of standard passbands, together 
with the size of our stellar sample, allowed us to derive the ba- 
sic physical properties of Fornax with high accuracy and mea- 
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Table 1. The journal of observations 



Nt. 


ir~\ 
ID 


Filter 


t-cxp 


A 


r YV 1 1 IV 1 J 


19 Oct. 


1995 


A 


B 


3x900 


1.02 


1.5 


19 Oct. 


1995 


A 


I 


3x900 


1.11 


1.5 


19 Oct. 


1995 


A 


V 


3x900 


1.06 


1.7 


20 Oct. 


1995 


B 


B 


3x1200 


1.09 


1.5 


20 Oct. 


1995 


B 


I 


1200 


1.01 


1.3 


19 Oct. 


1995 


B 


V 


600 


1.26 


1.4 


20 Oct. 


1995 


B 


V 


2x600 


1.22 


1.4 


21 Oct. 


1995 


C 


B 


3x1200 


1.25 


1.5 


21 Oct. 


1995 


C 


I 


3x1200 


1.03 


1.3 


20 Oct. 


1995 


C 


1 


240 


1.11 


1.0 


21 Oct. 


1995 


C 


V 


3x600 


1.10 


1.6 


20 Oct. 


1995 


D 


B 


1200 


1.34 


1.5 


21 Oct. 


1995 


D 


B 


2x1200 


1.03 


1.5 


21 Oct. 


1995 


D 


I 


3x1200 


1.01 


1.3 


21 Oct. 


1995 


D 


V 


3x600 


1.11 


1.6 


21 Oct. 


1995 


BKG 


B 


1800 


1.15 


1.6 


21 Oct. 


1995 


BKG 


I 


1500 


1.24 


1.4 


21 Oct. 


1995 


BKG 


v 


900 


1.36 


1.8 



sure details in its CMD that trace the less numerous populations 
and faster evolutionary phases. 

In particular, the B band turned out to be invaluable for 
studying the hot stars, be they old or young, whereas the wide 
baseline of the (B — I) color provides the best resolution of 
the different evolutionary phases in the color-magnitude dia- 
grams (cf. Smecker-Hane et al. 1994; Held et al. 



1999). Also, 



the availability of a comparison field allowed us to estimate the 
foreground and background contamination. The present pho- 
tometry will be the input database to model the star formation 
history (SFH) of Fornax (Held et al., in preparation) using pop- 
ulation synthesis techniques. 

The paper is organized as follows. The observations and 
data reduction are presented in Sect. ||. In Sect. ^] we present 
V, (B — I) color-magnitude diagrams of the Fornax field stars 
and discuss several interesting features with the help of theoret- 
ical isochrone fitting. The V, I luminosity function is derived 
in Sect. 



. 1 and used to estimate the distance to Fornax. This 
is confirmed by an independent distance estimate based on the 

we com- 



V luminosity of old-HB stars (Sect. 4.2). In Sect. 4.3 
pare the CMD of Fornax with template globular cluster RGB 
sequences using the standard (V — /) colors, and discuss the 
mean abundance and age of the dominant population. The color 
distribution of red giant stars in Fornax is analyzed in detail in 
Sect. %A_ , where the size of an intrinsic abundance spread is dis- 
cussed. Some light on the nature of the wide RGB of Fornax is 
shed by a co mpar ison of the spatial gradient of different popu- 
lations (Sect. 4.5). Our results and conclusions are summarized 
in Sect. || 

2. Observations and data reduction 



Fig. 1. The central area of Fornax reproduced from the Palo- 
mar Digital Sky Survey. The squares indicate the 10.'7 x 10.'7 
regions studied in this paper. The Fornax globular clusters #3 
and #4 are also indicated 

2.1. Observations 

The Fornax galaxy was observed on October 19-21, 1995 using 
DFOSC at the ESO/Danish 1.54m telescope. The detector was 
a 2048 x 2048 Loral CCD with pixel 0'.'40, covering a field 
of view of 13.'6 x 13f6. Due to non-uniform sensitivity near 
the edges of the CCD, the images were trimmed to a useable 
area of 1600 x 1600 pixels (i.e. 107 x 10.7). CCD readout 
by amplifier B in high-gain mode yielded a noise of 7.2 e~/px 
(rms) and a conversion factor of 1.31 e~/ADU. 

We observed 4 slightly overlapping fields in Fornax, plus 
one control field. A map of the location of the fields is shown in 
Fig. [I] The innermost field (region C) is centered on the glob- 
ular cluster #4, i.e. at about 4 arcmin from the galaxy center 
(as defined by SHS98). The outermost field (A) is located at 
~ 13'5 from the center. The journal of the observations is re- 
ported in Table |IJ The columns give the night, an image iden- 
tifier, the filter, exposure time and airmass, and the FWHM of 
the point spread function (PSF). The seeing was only fair, yet 
adequate to measure the relatively bright stars in our database. 
The third night had the most stable weather conditions. Several 
short exposure images, not included in this table, were used for 
checking the photometric zero points. 

2.2. Reduction and photometry 

The image processing was carried out with the ESO/MIDAS 
package in a standard wa y. Red uction follows the proc edures 
detailed by Saviane et al. ( 1996j Paper I) and Held et al. ( |l999 , 
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Paper II). For each field/filter combination, master images were 
produced by registering and coadding the long exposure im- 
ages. The PSF was not significantly degraded by this process. 
Stellar photometry was performed using DAOPHOT and ALL- 



STAR (Stetson 1987). The final PSF star catalogs contained 



~ 50 stars, and the best fit was obtained by fitting a Moffat 
(f3 — 1.5) function with a quadratic dependence on the x, y 
star coordinates. ALLSTAR was run twice on the sum images. 
In the second run, the star subtracted frames were searched for 
faint undetected objects that were added to the input lists of 
stars. The master photometric catalogs were created using the 
complete lists of stars as new inputs to ALLSTAR. 

2.3. Calibration 



Observations of Landolt's (1992) standard star fields were used 



to calibrate the photometry. The raw magnitudes were first nor- 
malized according to the following equation 



''ap 



2.5 logO 



cxp 



At) — k\ X 



(1) 



where m ap are the instrumental magnitudes measured in a cir- 
cular aperture of radius R = 6'.' 9, A t is the shutter delay and 
X the airmass. A shutter delay of —0.11 s was estimated from 
a sequence of images taken with increasing exposure times. 
The extinction coefficients ks = 0.235, ky — 0.135 and 
kj = 0.048 were adopted from the Geneva Observatory Pho- 
tometric Group data. The normal ized in strumental magnitudes 
were then compared to Landolt's ( 1992 ) values, and the follow- 
ing relations were found: 



B 
V 
V 
I 



= b' + 0.207 (B - V) + a B 
= v' + 0.0544 (B — V) + a v 
= v + 0.0489 (V - I) + c v 



0.00658 (V-I) 



ft/ 



(2) 
(3) 
(4) 
(5) 



where a B = 23.041, 23.025 and 23.038 for the nights 1, 2 
and 3, respectively. In the same order of nights, the other coef- 
ficients are a v = 23.774, 23.763 and 23.772; c v = 23.772, 
23.762 and 23.771; and finally a/ = 23.070, 23.054 and 
23.077. The standard deviations of the residuals were 0.018, 
0.013, and 0.022 mag in B, V (both equations), and I respec- 
tively. 

The PSF magnitudes were scaled to aperture magnitudes 
by assuming that m ap = mpsF + const. (Stetson 1987 ). Aper- 



ture magnitudes were measured for a sample of bright, isolated 
stars, and for each star we computed the difference with re- 
spect to the PSF magnitude measured on the coadded frames. 
The same reference aperture used for the standard stars was 
employed. The internal calibration uncertainty due to "aper- 
ture correction", estimated from the consistency of the zero 
point derived from several individual images, is of the order 
0.01 mag for all filters. 

The instrumental magnitudes and colors for the stars ob- 
served in at least two filters were calibrated either with an it- 



erative procedure or by solving a system of 3 equations of the 
form 

TO st d = m in st + k m color std + a m 

where m s td and colored are the magnitude and color in the 
standard system, and mj nst are the instrumental magnitudes. 

After independent calibration, we performed a verification 
of the photometric zero points of the catalogs in the 4 zones by 
comparing stars in the overlap strips. Since the mean system- 
atic deviations between BVI magnitudes measured in the field 
C and D (both observed during the third, most stable night) 
are less than 0.03 mag, we chose to refer all photometry to the 
zero point of the central field C, and applied small zero-point 
corrections to our photometry in the field A, B, and D. We con- 
servatively adopt an uncertainty of 0.03 mag as our systematic 
error in all bands. 



2.4. Comparison with previous studies 

As a further check of the accuracy of our photometric zero 
point, we compared our results with previous data in the liter- 
ature. The onl y pub lished photometry tables are those of Buo- 
nanno et al. ( 1985 ). Their Tables 6 and 7 report the values 
of V and (B — V) for all the stars measured in two separate 
2x2 arcmin 2 fields, called Al and A2, which are included 
in our fields C and A. The two sets of measurements for the 
A2-A field pair are in good agreement. The median differences 
(this paper - Buonanno et al.) are —0.013 in V and —0.017 in 
B — V, with standard deviations of 0. 16 mag in both cases. 

The consistency between the zero-points for the Al-C pair 
is still good, yielding median residuals 0.014 (er =0.20) in V 
and 0.024 (0.25) in B - V. 



2.5. Artificial star tests 

Extensive artificial star simulations were performed to evalu- 
ate the uncertainties of our photometry and the completeness 
of the data. The simulations were carried out for the fields A 
and C, which represent the lowest and highest crowding in our 
frames. A list of input stars was created for each V master im- 
age, with uniformly distributed magnitudes. The star coordi- 
nates were generated over a grid of triangles with a small ran- 
dom offset from the vertices, a configuration allowing to add 
the largest number of non-overlapping simulated stars. Error 
estimates based on randomly placed artificial stars may not be 
realistic if there is a significant amount of clustering among 
real stars. This caveat does not seem to apply to our relatively 
uniform Fornax fields, though. 

The same artificial stars were used in all bands, using ran- 
dom (B — V) and (V — 7) colors so that the stars were uni- 
formly distributed in the color-magnitude diagrams. We typ- 
ically added ~ 36000 stars per filter in 80 experiments. The 
frames with the artificial stars were then reduced using exactly 
the same procedures as for the original images. For each fil- 
ter, the retrieved artificial stars were matched to the input list 
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Fig. 2. The color-magnitude diagram 
of Fornax in the (B — I), V plane. This 
diagram includes about 42500 stars in 
all fields. The most noticeable fea- 
tures are the wide RGB made of old 
and intermediate-age stars, the upper 
AGB tail, a young main sequence, and 
a prominent red clump together with 
fainter, older HB stars. The blue main 
sequence clearly merges into a mix 
of subgiant branches. The outlined re- 
gions have been used for counting stars 
in different evolutionary phases (see 
text) 



Table 2. The photometric errors from artificial star experi- 
ments 



V,B,I 


<j\- 




cri 


(j v 




<7l 






A 






c 




14.75 






0.003 








15.25 






0.003 








15.75 






0.003 






0.011 


16.25 






0.004 






0.012 


16.75 






0.005 






0.014 


17.25 




0.013 


0.006 




0.010 


0.017 


17.75 




0.014 


0.013 


0.014 


0.014 


0.018 


18.25 


0.013 


0.014 


0.017 


0.016 


0.019 


0.019 


18.75 


0.017 


0.019 


0.020 


0.020 


0.018 


0.026 


19.25 


0.021 


0.024 


0.028 


0.025 


0.023 


0.032 


19.75 


0.024 


0.031 


0.037 


0.032 


0.030 


0.045 


20.25 


0.030 


0.038 


0.049 


0.040 


0.036 


0.064 


20.75 


0.038 


0.052 


0.069 


0.061 


0.049 


0.082 


21.25 


0.057 


0.069 


0.098 


0.081 


0.064 


0.108 


21.75 


0.073 


0.085 


0.144 


0.103 


0.082 


0.156 


22.25 


0.107 


0.105 


0.150 


0.136 


0.106 




22.75 


0.116 


0.134 


0.181 




0.124 





by means of their coordinates. The stars recovered in differ- 
ent colors were matched, and the raw catalog calibrated just as 
the original photometry. The standard deviations of the mea- 
surement errors Am, calculated in 0.5 mag bins, are given in 
Table |2|. The first column gives the bin centers, columns 2 to 
4 and 5 to 7 list the errors obtained for the fields A and C, 
respectively. The measured standard errors span a range from 



~ 0.01 mag near the tip of the RGB to ~ 0.15 mag close to 
the limiting magnitudes. Errors and the limiting magnitudes are 
consistent with the different crowding conditions and exposure 
times in the two fields. Contour plots of the completeness lev- 
els were produced by dividing each CMD in cells with color and 
magnitude steps of 0.5 and 0.2 mag, respectively, and counting 
the artificial stars within each cell before and after reduction. 
We used the same acceptance criteria as for the galaxy's CMD, 
i.e. a star was counted in a cell only if it was recovered in each 
of the 3 filters. The completeness array was computed as the 
ratio between the post-reduction and the input simulated star 
counts, and median filtered with a 5 x 5 box to reduce the noise 
in the contour plots. Examples of the resultant completeness 
contours will be shown in Sect. ||. 

3. Color-magnitude diagrams 

The combination of a wide photometric baseline and large sam- 
ple size employed in this study allowed us a very detailed view 
of the evolved stellar populations in Fornax. Fig. || presents 
the V, (B — I) color-magnitude diagram of this nearby dwarf 
spheroidal, showing an excellent separation of stars in different 
evolutionary phases. We now briefly describe the many inter- 
esting features seen in this CMD. 

- A wide red giants branch, comprising stars older than ~ 1 
Gyr. The color spread is much larger than expected from 
photometric errors (cf. Sect. p~4| ). While the red side of the 
RGB shows a sharp edge, the stars are spread on the blue 
side forming sparsely populated sequences distinct from 
the RGB (cf. SHS98). This is more evident in Fig. (| Note 
that the foreground and background contamination is virtu- 
ally negligible in the relevant regions of the diagram. The 
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Fig. 3. An enlarged view of the color- 
magnitude diagram of Fornax showing 
the features produced by core helium- 
burning stars, and morphological de- 
tails of the red giant branch. Two 
isochrones from the post-MS models of 
Bertelli et al. (11994 with Z=0.004 and 



ages 300 and 400 Myr (top to bottom) 
have been superimposed to the data. 
The plume of stars above the red clump 
is composed of intermediate mass stars 
burning helium in the core after leaving 
the young main sequence. 



color-magnitude diagram of foreground and background 
objects in our control field, covering an area 10'.7 x 10.'7, 
contains only 154 objects down to V ks 22.5. 
Above the RGB tip at V ~ 18.4, there is a well-developed 
upper AGB tail extending to colors as red as (B — V) and 
(V - I) ~ 3, or [B - I) ~ 6. The upper AGB consists 
of intermediate-age C and M stars, the latter comprising a 
small group of stars just above the RGB tip (see SHS98). 
A rich red clump contains the majority of HB stars of a nu- 
merous intermediate-age and metal-enriched population; in 
the following we will re fer to it simply as red clump (RC) 
(cf. Demers et al. 



1994 



SHS98). About 0.2 mag fainter, 
a horizontal branch originating from an older population 
is clearly seen, indicated in the following as "old HB" (or 
HBold)- We also notice the instability strip mostly popu- 
lated by RR Lyrae variables, whose random phase colors 
and magnitudes define a band Sa 1 mag thick. RR Lyrae 
variables are present in the Fornax field (Buonanno et al. 



1985; SHS98; Mateo 1998). Blue HB stars are hardly seen 



in this diagram. If any exist, they are confused with the 
young main sequence stars. 

A blue plume reaching V ~ 20, identified with a young 
main sequence of ~ 0.1 Gyr old stars (Beauchamp et al. 
1995| ; SHS98). 

An almost vertical plume originating from (B — I) ~ 1.8, 
V ~ 21, i.e. just above the red clump, extending up to 



V ~ 19.3. They are, as it will be shown below, core helium- 
burning stars with mass in the range ~ 2 Mq. 

- At fainter magnitudes, the main sequence merges into a 
heavily populated region ~ 1 mag below the HB, involv- 
ing a mix of subgiant branches of different ages. 

- We also note the small clump of stars at V ~ 20.4 on the 
red giant branch, an example of the short-lived evolutionary 
phases that can be revealed in stellar populations using ad- 
equately large data samples. This feature is identified with 
the AGB bump, a clumping of stars due to a slowing down 
of the luminosity incre ase at the beginning of AG B evolu- 

Alves & Sarajedini |1999| ). We will 



1998 



tion (e.g., Gallart 
return on this point in Sect, kkl 



Overall, this diagram shows that Fornax went on forming stars 
from an early epoch (> 10 Gyr ago) until about 100 Myr 
ago. A closer picture of the CMD is presented in Fig. |L show- 
ing details of the core helium-burning stars, which are im- 
portant diagnostics of stars formation histories in galaxies. 
The prominent red clump comprises stars with age 2-10 Gyr 
and mass approximately 0.9 to 1.4 Mq. Its mean luminosity 
and color bears information on the mean age of Fornax. We 
will return on this point later. A comparison with theoretical 



isochrones of Bertelli et al. (1994) allows one to establish the 



nature of the stars producing the plume above the RC. These 
are intermediate-mass core helium-burning stars (2.4-2.9 Mq), 
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counterparts of the young main sequence stars in the age range 
0.3-0.5 Gyr, which started burning helium in a non-degenerate 
core. 

This stage (also known as the "blue-loops") represents an 
important indicator of the metallicity of the young population 
and therefore of the che mical enrichment histor y of galax- 
ies (e.g., Aparicio et al. 1996 ; Cole et al 1999| ; and refer- 
ences therein). Fig. ^] shows that fitting the RC plume requires 
isochrones having [Fe/H]^ —0.7, i.e. significantly more metal- 
rich than the bulk of the Fornax stellar population. This rel- 
atively high metallicity of the younger stars may probably ex- 
plain why Fornax seemingly lacks a large population of anoma- 
lous Cepheids ( AC), which are so numerous in the metal-poor 
dS ph Le o I (see discussion of the instability strip in Caputo et 
al. |1999 ). Searches of AC's in Fornax are underway (Bersier & 
Wood|l999|). 



The color-magnitude diagrams for the innermost and outer- 
most region in this study, shown in Fig. ^j, illustrate the remark- 
able variation of the stellar populations in Fornax with galacto- 
centric radius. The inner field (C, bottom panel) shows all of the 
features noticed in the total CMD. The young main-sequence 
stars are less numerous in the outer field, even accounting for 
the lower stellar surface density, and there are very few blue 
stars brighter than V ~ 22. The young core-He-burning stars 
(the plume above the RC) follow the trend of the blue main se- 
quence stars. In contrast, it is interesting to note a hint of a blue 
horizontal branch in the outer field. Together with the detec- 
tion of RR Lyrae stars, our data provide evidence for a small 
old, metal-poor field population similar to that of the Fornax 
globular clusters. A population II halo seems to be common not 
only in dwarf spheroidals (see Mate o 1998 ) but also in dw arf 
irregulars (e.g., Minniti et al. 1999 ; Aparicio et al. 1997 ). A 
quantitative estimate of the population gradient in Fornax will 
be given in Sect. 4.5. 



4. Analysis and discussion 

4.1. Luminosity function and distance 

The red giant luminosity function (LF) and distance to Fornax 
was derived using star s within ±2cr from the fiducial sequence. 
As shown in Sect. 4.4 , this implies selecting the dominant stel- 
lar population in Fornax. Luminosity distributions were ob- 
tained both in V and in I by counting stars in 0.2 mag bins 
down to below the red clump. Since at these bright magni- 
tudes our photometry is virtually complete, there was no need 
to correct the observed LF's for incompleteness. Foreground 
and background contamination is not a concern, either, because 
the number of field objects in the proximity of the RGB is neg- 
ligible. 

The cutoff in the / luminosity function corresponding to 
the maximum luminosity reached by red giants before they ig- 
nite the He burning, has proven to be a good distance indicator 
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(Madore & Freedman 1995; see also Salaris & Cassisi 1998). 
We measured the / magnitude of the RGB tip separat ely fo r 



our 4 Fornax fields, following the methods of Lee et al. ( Jl993[ ) 



Fig. 4. A comparison of the color-magnitude diagrams in an 
outer and inner region of Fornax. Field A (top panel) is located 
about 14' from the center of Fornax, while field C (bottom) 
samples an area near the galaxy center. The contour lines rep- 
resent the 30%, 50%, 70% and 90% completeness levels of our 
photometry. Note that the completeness in the crowded field C 
is lower than in field A. The number of young main sequence 
stars and their helium-burning counterparts decreases from the 
inner to the outer field, allowing perceiving the bluer HB stars 
belonging to the oldest Fornax field population 



The mean of the four values is < Itrgb >= 16.72 ± 0.10. 
Although the scatter of the individual measurements is small 
(~ 0.02 mag), we have adopted a larger uncertainty to take into 
account both the intrinsic precision of the tip detection method, 
which is about a half of the 0.2 mag bin, and the zero point un- 
certainties. 

We then calculated the distance to Fornax using the rela- 
tions of Da Costa & Armandroff ( 1990| ). This procedure, often 
applied to derive the distance of dwarf galaxies, implicitly as- 
sumes that the / magnitude of the tip is little affected by age. 
Theoretical models indeed show a dependence of the tip lumi- 
nosity on the age, but this is more pronounced for very metal- 
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models (Bertelli et al. 1994 ). We have measured the location of 



Fig. 5. The V RGB luminosity function of Fornax, showing 
the sharp cutoff at the RGB tip, the prominent red clump at 
My ss +0.4, and the "AGB bump", a clumping of intermediate 
age stars at the beginning of their AGB double-shell burning 
phase (M v m -0.4) 



the AGB bump in Fornax by performing a Gaussian fit to the 
LF in the region of the bump. We found V = 20.40 ± 0.04, 
where the error is mainly set by the zero point uncertainty. 
This corresponds to a luminosity My ~ —0.39 ± 0.04, with 
an additional 0. 1 mag uncertainty on the Ay extinction. Sim- 
ilarly, we measured a mean V magnitude for the red clump 
Vrc = 21.18 ± 0.04 mag, corresponding to M v = +0.39. 
Thus the detected clump is 0.78 ± 0.06 mag brighter in V than 
the red HB clump. For an assumed age of 5 Gyr and the mean 
metallicity of Fornax, these measurements confirm the identi- 
fication with the AGB bump and rule out alternative identifica- 
tions with the RGB bump. The RGB bump is expected to be 
near the HB level fo r a me tallicity [Fe/H] w — 1 and age 5 Gyr 
(Alves & Sarajedini 1999 ). In Lyndsay 113, a 5 Gyr old clus- 
ter in the Small Magellanic Cloud having metallicity compara- 
ble with that of Fornax, the RGB bump is found to be ~ 0.15 
mag b righter than HB stars in the same cluster (Mighell et al. 
1998). Clearly, stars in the RGB bump will be outnumbered 



by the overwhelming red clump. Observational data like those 
presented here for the Fornax dwarf are important to constrain 
evolutionary models, which in turn are necessary to interpret 
the stellar population of Local Group galaxies. 



poor populations ([Fe/H]< 
(e.g., Caputo et al. |1999|). 



1.7) and ages younger than 5 Gyr ^ 2 Distance based on the old horizontal branch 



The relations of Da Costa & Armandroff ( Jl 99Cj > give the / 
bolometric correction as a function of color of the stars near the 
RGB tip, and the bolometric luminosity of the tip as a function 
of metallicity. The dereddened color of the RGB tip in For- 
nax, calculated as the median (V — /) within 0.1 mag from 
the tip, is (V — -Oo.trgb = 1-59 ± 0.06, where the error re- 
flects the scatter of the values obtained in our four fields plus 
the absolute zero point uncertainty. We adopt ed a reddening 
E B -v = 0.03 ± 0.03 from Burstein & Heiles ( |1982[ ), yielding 
£V-j = 0.038 + 0.038 and Ai = 0.058 + 0.058. The bolomet- 
ric correction is then BCi = 0.495 ± 0.015, while a metallic- 
ity [Fc/H] = -1.39 ± 0.15 (cf. Sect. O) implies M™ GB = 
-3.55 ± 0.01. We thus obtain Mj RG ®~= -4.04 ± 0.02, and 
a distance modulus (m — M)o = 20.70 ± 0.12, corresponding 
to 138 ± 8 kpc. 

Previous distance estimates range from (m — M)o = 
20. 59 ± 0.22 (Buonanno et al |1985| ) to 20.76 (Dem ers et 
al. |1990| ; Buonanno et al. |1999| ). Sagar et al. ( |1990| ) found 
(m— M)o = 20.7. The present estimate therefore confirms ear- 
lier results. This value is also consistent with th e dista nce mod- 
uli of Fornax globular clusters (Buonanno et al. |1998[ ), yielding 
an average modulus (m — M)q ~ 20.62 ± 0.08. 

Using this distance estimate, we plot in Fig. |]the V lumi- 
nosity function of the red giant stars in the inner region of For- 
nax. Besides the obvious red clump at My ~ 0.4, we notice the 
small yet significant peak near V — 20 .4 that we identify with 
the AGB bump, a clustering of stars that occurs at the beginning 
of helium shell-burning evolution. Gallart ( 1998| ) has recently 
discussed the presence of this feature in the LMC and M31 
where its location agrees with the prediction of stellar evolution 



An independent estimate of the distance to Fornax was ob- 
tained from the mean level of its old- HB field stars. The mean 
magnitude of the HB was measured by fitting a Gaussian to the 
V mag distribution of the stars in the range 21.2 < V < 21.7, 
1.1 < B — I < 1.5. We necessarily included only the red part 
of the HBold, since the bluer horizontal-branch stars appear 
to be mixed with the blue stars on the young main sequence. 
Note, however, that we do not include any RC stars (which are 
brighter than the HBold and RR Lyrae variables). The mean 
level of the red HB is V H b = 21.37 ± 0.04, where the uncer- 
tainty reflects the scatter of the values measured in the different 
fields (larger than the formal error on the mean magnitude), and 
the systematic error of the V zero point. Buonanno et al. ( 1998) ) 
found a mean level Vhb = 21.25 ± 0.05 for the HBoLDof four 
globular clusters in Fornax. 

Using this value for V H b, and Ay = 3.2 E(B — V) = 
0.096, we calculated the distance modulus of Fornax on the Lee 



et al. (1990) distance scale, using their relation for the absolute 
visual magnitude of RR Lyrae variables, 



My = 0.17 [Fe/H] +0.82 



(6) 



for a helium abundance of Y = 0.23. 

Assuming for red HB stars the relatively metal-rich nom- 
inal metal content of RGB stars, [Fe/H]« —1.4, this rela- 
tion would give My K = 0.59 mag and a distance modulus 
(m — M)o = 20.69 ± 0.04 for a population with age com- 
parable to that of Galactic globular clusters. This uncertainty 
includes internal and photometric errors only. However, the 
mean metallicity of the old-HB stars is probably lower. If the 
old population in Fornax is relatively metal-poor, of the order 
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(V-l) 

Fig. 6. A comparison of our total color-magnitude diagram of 
Fornax with the giant branches of templ ate G alactic globular 
clusters from Da Costa & Armandroff (1990), scaled to the 



distance and reddening of Fornax. The globular clusters span 
a metallicity range from [Fe/H] = -2.2 to [Fe/H] = -0.7 



[Fe/H]« —1.8 ( as we suggest in Sect. 4.5), the relation given by 
Lee et al. ( |1990| ) would imply M$ R = 0.51 mag and a distance 
modulus (to - M) = 20.76 ± 0.04. The level of the red HB 
(distinct from the clump) is probably the result of contributions 
from stars in a range of ages and metallicities. For this reason 
we refrained from applying any uncertain correction to convert 
the measured mean magnitude of the red HBold to an equiv- 
alent magnitude of RR Lyrae variables. Further, this distance 
modulus based on the HB level is affected by the uncertainties 
on the luminosity of HB stars as a function of age and metal- 
licity. A discussion of the alternative distance scales, however, 
is beyond the scope of this paper. 

This measurement of the distance to Fornax based on its old 
horizontal branch star luminosity confirms the distance modu- 
lus estimated from the RGB tip. This consistency is not unex- 
pected, since both the RGB tip method of DA90 and the HB 
absolute magnitude obtained for the HBold are based on the 



distance scale of Lee et al. (199C). These two distance mea 



surements use / and V magnitudes, respectively, which are ob- 
servationally independent. 

4.3. Mean abundance and age 

The mean metal abundance of the bulk of the Fornax popula- 
tion was derived by direct comparison of the red giant branch 
in the /, (V — I) color-magnitude diagram with the ridge lines 
of globular clusters from Da Costa & Armandroff ( 1990 ) (see 
Fig. ||). Our procedure is fully described in Paper I and II, and 



is only briefly outlined here. In short, we calculated the aver- 
age color shift, 6(V — I)o, between the Fornax RGB and the 
Galactic cluster fiducial loci. An interpolation of the relation 
between the mean color shifts and the globular cluster metal- 
licities (actually a quadratic fit) provides an estimate of [Fe/H] 
for the dwarf spheroidal. This procedure was applied to the 2a- 



selected RGB sample (cf.Sect. 4.1), in two luminosity intervals 
(-4.0 < Mi < -3.0 and -3.0 < Mj < -2.0), yielding a 
metallicity [Fe/H]= -1.45 ± 0.11 and [Fe/H]= -1.33 ± 0.15 
dex, respectively. The mean of the abundances determined in 
these two magnitude bins was adopted as our final estimate. 
The resultant value, [Fe/H]= —1.39 ± 0.15, is in good agree- 
ment with previous work. We find no evidence for a metallicity 
gradient among the regions studied here, to within the errors. 

However, the measurements of mean abundance based on 
the color of the RGB are subject to the well-known difficulty in 
disentangling the effects of age and metallicity on the effective 
temperature of red giant stars (the "age-metallicity degener- 
acy"). Thus we need to estimate the mean age of Fornax before 
discussing further its mean metal abundance. When compared 
with the predictions of stellar evolution models (e.g., Bertelli et 
al. 1994; Caputo et al. 1995), the position of core He burning 



stars in color-magnitude diagrams may provide a useful age in- 
dicator (e.g., Caputo et al. 1999; Girardi 1999; and references 
therein). The RC comprises core helium-burning stars of dif- 
ferent ages, so that its location bears information on the mean 
age of the intermediate age population, weighted by the age 
distribution function. Thus, similarly to what we had done for 
the HB, we measured the mean (B — I) color in addition to the 
V luminosity for the red clump. The mean magnitude, already 
reported above, is Vac = 21.18 ± 0.04 mag, corresponding 
to My — +0.39, in excellent agreement with Demers et al. 



1994). This means that the RC is 0.19 ± 0.06 mag more lu- 



1998), but less 



1999a). The mean color 



minous in V than the old HB stars, a value that appears con- 
sistent with the difference in age of a 13 Gyr old popu lation 
and a 5 Gyr old bulk component (see Caputo et al. 1999 ). The 
clump is quite extended in luminosity (~ 0. 6 mag) , compa- 
rable with that of Carina (Hurley-Ke ller et al. 
than that of Leo I (cf. Gallart et al. 
is < B — I >rc= 1-79 ± 0.04. The uncertainties include 
the field-to-field scatter, comparable with the photometric mea- 
surement errors, and the zero-point uncertainty. The relation 
(V - 1) = 0.457 (B -I) + 0.147, obtained from a linear fit to 
the color-color relations for the Fornax red giants in the range 
1.0 < (B - I) < 3.5, yields < V - I > RC = 0.965. This 
value shows excellent agreement with the results of Buonanno 
et al. ( 1999| ). By fitting a parabola to the fiducial points of the 
RGB, we estimated the interpolated RGB color at the RC level 
(V— J ~ 1.07 mag), a value also confirmed by inspection of the 



WFPC2 color-magnitude diagram (Buonanno et al. 1999). The 



difference in color between the red clump and red giant stars 
at the same luminosity is then <W-j),bc = 0.10 mag, with 
an estimated uncertainty of 0.02 mag. This result can be com- 



pared with the model predictions of Girardi (199S; and priv. 
comm.) based on the models of Girardi et al. ( 1999 ), which 
are in accord with the empirical calibration of Hatzidimitriou 
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Table 3. Observed and instrumental color dispersions am-i) 
for the main component of the RGB color distributions in fields 
A and C. 



I 






tJA(instr) 


(Tc(instr) 


O"0,A 


O"0,C 


17.20 


0.082 


0.067 


0.024 


0.029 


0.078 


0.060 


18.20 


0.063 


0.077 


0.036 


0.042 


0.052 


0.065 


19.20 


0.088 


0.093 


0.061 


0.073 


0.063 


0.058 



( |1991| ). For a metallicity Z=0.001 (but <5(v_/) iRC is relatively 
independent of abundance for metal-poor populations) our re- 
sult is consistent with a mean age of the order 5. 4 ± 1 .7 Gyr. 
This value is close to the estimate of Sagar et al. ( |1990| ), based 
on best fitting of Yale isochrone s, and definitely larger than the 
age estimated by Demers et al (1994). Most interestingly, the 
mean age obtained from the clump location appears to be con- 
sistent with the presence of MS evolved stars in the same age 
interval, as observed with HST (Buonanno et al. 1999). This 



results is quite encouraging for application of this age indicator 
to more distant Local Group galaxies, whose main-sequence 
turnoff cannot be directly measured. 

If we now assume a mean age of approximately 5 Gyr for 
the bulk of the Fornax stars, the observed RGB color would 
imply a metallicity significantly larger than the formal result 
obtained above from a comparison with globular clusters. We 
have estimated the effects of age by comparing theoretical 
isochrones of different ages (e.g., 5 and 15 Gyr) at a given 
metallicity (from Bertelli et al. 1994). By measuring the (V — I) 
colors at Mj = —2.5 predicted by model isochrones with 
Z=0.001 ([Fe/H]= -1.3), we find that a 5 Gyr isochrone is 
bluer by ~ 0.09 mag than a 15 Gyr model isochrone. This 
effect mimics a metallicity difference of ~ 0.4 dex using the 
fiducial loci of globular clusters (cf. Paper II; Caputo et al. 
1999| ; Gallart et al. |1999a| ). Thus, if the body of Fornax stars 
is ~ 5 Gyr old, the measured location of the peak of the RGB 
is necessarily indicative of a higher mean metallicity, of the or- 
der [Fe/H]= —1.0 (clearly the correction is somewhat model 
dependent). We regard this value as the most appropriate esti- 
mate of the mean metal abundance of the dominant stellar pop- 
ulation in Fornax. With this correction, the Fornax metallicity 
turns out to be very close to that of Sagittarius, a dSph whic h 
has a comparable total luminosity (e.g., Bellazzini et al. 1999 ). 



4.4. The color distribution of Fornax red giants 

One of the main results of this paper, made possible by the 
size of our stellar sample and photometric baseline, is a de- 
tailed analysis of the color distribution function ( CDF) of the 
red giant stars in Fornax. Fig. [7] shows the distribution of the 
(B — I) color residuals about a preliminary fiducial sequence, 
in the magnitude range 17.7 < I < 18.7 (-3 < A// < -2), 
for the inner and outer field. While these histograms confirm 




Fig. 7. The color distribution of the red giant stars in Fornax, 
plotted separately for the inner (bottom panel) and outer region 
(top panel). The histograms represent the distributions of the 
color residuals of individual stars from a median RGB fiducial 
sequence, in the magnitude interval 17.7 < I < 18.7. Error 
bars represent Poisson errors. The color distribution is quite 
well fitted by the sum (continuous line) of two Gaussian func- 
tions (dashed lines) suggesting a two-component model for the 
metallicity (or age) distribution of Fornax stars. 
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Fig. 8. Plot of the measured color scatter am-i) of the princi- 
pal RGB component in our fields A and C (circles), compared 
with instrumental dispersions (triangles). 



the well-known wide color range of the RGB stars in Fornax 
(e.g., Buonanno et al. 



1985 Sagar et al. 199C; Grebel et al. 



1994; Beauchamp et al. 1995), they also show that the color 



distributions cannot merely be described using a single "color 
dispersion". Rather, the CDF is more appropriately described 
as roughly bimodal, showing a principal peak and a bluer com- 
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ponent extending to A(B — I) ~ —0.4. This color distribu- 
tion function is quite well modeled by the sum of two Gaus- 
sians. The main component of the distribution represents the 
bulk of the red giant population, a mix of old and (mostly) 
intermediate-age stars. The secondary peak is centered at about 
A B _/ = -0.20. 

On the other hand, we notice a relatively well-defined cut- 
off on the red side of the RGB, indicating the lack of any sig- 
nificant metal-rich component similar to the stellar population 
of 47 Tuc, or even less metal-rich if we assume a mean age 
younger than that of Milky Way globular clusters. This absence 
sets an important constraint for modeling the chemical enrich- 
ment of the Fornax dwarf. Both components are wider than ac- 
counted for by instrumental errors. The dispersions of the two 
components, in the luminosity range —3 < Mi < —2, are 
(7%_j — 0.063 (central peak) and cr B _ I = 0.088 (blue compo- 
nent) in the inner region, and <j b _j = 0.077, <J h B _j = 0.099 
in the outer field. Table || gives the observed dispersions for 
the main component of the Fornax field population for the in- 
ner and outer region in 3 magnitude intervals. Also given in 
Table | are the instrumental errors a(B — I) obtained by fit- 
ting a Gaussian to the color residuals of artificial stars, ex- 
actly in the same way as for the real data. The observed and 
instrumental dispersions are also compared in Fig. |J The in- 
trinsic (B — I) color dispersions, calculated as the quadratic 
difference between the observed and the instrumental scatter, 
are given in the last two columns of Table ||. In the luminos- 
ity interval — 3 < Mj < — 2 the intrinsic color spread of the 
main RGB population is <jq(B — I) = 0.06 ± 0.01 mag. Us- 
ing again the color-color relations for the Fornax red giants to 
convert (B — I) color spreads into equivalent dispersions in 
(V — I), and the calibration of RGB color shifts as a function 
of metal abundance variations, we obtained a metallicity spread 
C[Fe/H] = 0.12 ± 0.02 dex for the dominant field population. 

Then the (2a) metallicity range for the bulk population of 
Fornax would be approximately —1.65 < [Fe/H] < —1.15, or 
— 1.25 < [Fc/H] < —0.75 if a correction for the mean age is 
applied. This intrinsic metallicity range is significantly lower 
than the abundance spread quoted by most previous studies for 
the red g iant branch as a whole. According to Beauchamp et 
al. (1995), the total range i n [Fe/ H] is 0.8 dex, com parabl e to 
that found by Sagar et al. ( |1990| ) and Grebel et a l. (]1994| ), A 
smaller spread (^ 0.1 dex) was found by Geisler (1994). This 
discrepancy probably results from the coarser metallicity reso- 
lution of the colors employed in past studies, with the notable 
exception of the Washington colors of Geisler (1994). We con- 
clude that a small abundance spread seems in fact more ap- 
propriate to describe the intermediate-age field population in 
Fornax. While the metallicity dispersion given above is com- 
parable to that of Leo I (e.g., Gallart et al. 1999a) ), it appears 
to be smaller than the abundance spread found in the majority 
of dwarf spheroidal galaxies (Da Costa |1998| Mateo |1998[). In 



et al. 1999). Had we considered our Fornax RGB color distri- 



a few cases, wide range in metallicity has been confirmed by 



low- and high- dispe rsion spectroscopy (e.g., Cote et al. 1999; 
Shetrone et al. 1998 ). The mean value of the metallicity spread 
for Galactic dSph and satellites of M 3 1 is 0.37±0.03 dex (Cote 



bution as a whole, we would have obtained a metallicity spread 
of the same order (-2.0 < [Fc/H] < -0.7, ±2a range), in 
good agreement with previous studies. 

Since age is known to affect the RGB color, [Fe/H] dis- 
persions derived by the width of the giant branch should be 
taken with caution in view of a possible contribution of an 
age spread. As argued above, an age range of the order 5 Gyr 
(which is that of stars making up the main RGB) is sufficient to 
mimic a metallicity range of ±0.1 dex. Thus we might assume 
that the width of the RGB main component is entirely due to 
the age spread of its populations. The situation is more com- 
plex, though, and the effects of a metallicity and age range on 
the color distribution depend on the details of the star forma- 
tion and chemical enrichment history. Successive stellar gener- 
ations are expected to be progressively more metal-enriched, so 
that younger stars (implying a bluer RGB) will generally have 
higher metal abundance (leading to a redder RGB). The two 
effects - of a younger age and higher metallicity - will act in 
opposite directions, and may even compensate each other as it 
appears to be the case for Carina (Smecker-Hane et al. 1994). 



Similarly to Carina, the abundance spread we have found for 
the dominant population of Fornax may represent a lower limit 



(see also Paper II; Gallart et al. 1999a; for similar considera 



tions for other dSph's). This issue shall be more quantitatively 
investigated in a following paper. 

We return now to discuss the nature of the population mak- 
ing up the blue tail of the CDF, which is until now far from 
established. Qualitative examination of the CMD's is not suffi- 
cient to establish whether the blue tail of the CDF represents an 
old, metal-poor population, or is made up of young red AGB 
stars. However, we will show in Sect. 4.5 that there is defi- 



nite evidence that the bluer RGB stars are old and metal-poor, 
which implies that the extended color distribution shown in 
Fig. ^ can be interpreted as a metallicity distribution. In con- 
clusion, a model involving two populations seems to provide a 
good description of the star content of the Fornax dSph, with 
the older population having [Fe/H]= —1.82 with a dispersion 
of 0.20 dex, and the dominant, intermediate-age population 
with [Fe/H]w —1.0 ± 0.15. Our large-field data confidently 
rule out the presence of a distinct metal-rich population with 
abundance comparable to that of 47 Tuc, even accounting for a 
mean age of 5 Gyr for the Fornax bulk population. 

This complexity is common to most of the other dwarf 
spheroidals. For example, two distinct star formation epochs 
have been recently revealed in Sculptor by Majewski et al. 
(1999). In this galaxy, a detailed analysis of the RGB morphol- 



ogy showed the presence of two distinct RGB bumps consistent 
with the presence of a metal-poor population of [Fe/H]~ —2.3, 
and a population of [Fe/H]~ —1.5. Also the recent study of the 



star formation history of Leo I by Gallart et al. (1999b ) indicate 



that most of the star formation activity (80%) occurred between 
7 and 1 Gyr (mean 4 Gyr) while the contribution of the older 
phase was small. A wide metallicity range and a composite 
population, although with a higher mean abundance, has also 
been inferred in the Sagittarius dSph, a galaxy similar in many 
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Fig. 9. Radial trends in the fraction of young main sequence 
stars (filled circles), HB stars (triangles), stars on the blue-RGB 
(squares) and RGB (open circles) relative to the number of red 
clump stars. The logarithm of the ratios is plotted against the 
effective distance from the Fornax center. 

Table 4. Star counts normalized to 100 stars in the red clump 
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respects to Fornax (Bellazzini et al. |1999j). Also, the metallicity 
distribution of stars in the small elliptical M 32 shows a metal- 
rich peak ([Fe/H]~ —0.2) with a low-m etallici ty tail extending 
to about [Fe/H]'- -1.5 (Grillmair et al. |1996| ). It is also inter- 
esting to note the analogy with the extremely broad metallicity 
range fo und in the halo of the nearby elliptical NGC 5 128 (Har- 
ris et al. 1998 ), where the shape of the metallicity distribution 
suggested a two-phase in situ model. 



4.5. Population gradients 

A comparison of the color-magnitude diagram in the different 
regions in this study provided important clues regarding the 
origin of the stellar populations in Fornax, and in particular on 
the nature of its complex red giant branch. Were the bluer RGB 
stars old and metal-poor, one would expect a higher fraction of 
them in the outer fields, on the basis of the population gradient 
detected by SHS98. Clearly the opposite finding, i.e. a larger 
RGB blue tail in the inner regions, would indicate a connection 
to the more recent bursts of star formation. 

In order to measure the radial gradient in the stellar pop- 
ulations in Fornax, stars in different evolutionary phases were 
counted separately in different radial bins. The CMD regions 
chosen for counts include the red clump, the red part of the 
HBold, the blue-loop helium-burning stars (BL), the young 
main sequence (yMS), and the red giants (those in the main- 
stream giant branch, RGB, and in the bluer component, 6RGB). 

The reader is referred to the boxes outlined in Fig. ||. The 
results of star counts are presented in Table [I[ where we list the 
effective galactocentric distance and the percentage of stars in 
all the CMD regions relative to the number of RC stars. 

The fraction of young main sequence, old HB, blue-RGB 
and mainstream RGB stars are also plotted on a logarithmic 
scale in Fig. [| As previously noticed by SHS98, the young 
stars are more centrally concentrated than the dominant inter- 
mediate age component, indicating that recent star formation 
took place preferentially in the central regions. The counts on 
the RGB as expected follow those of RC stars. Conversely, the 
HB stars are preferentially found in the outer regions. 

Most importantly, the bluer RGB stars closely follow the 
radial trend of the horizontal-branch stars (Fig. ^|). This result 
unambiguously demonstrates that the sparse sequence on the 
blue side of the Fornax RGB belongs to the old and metal-poor 
population (^ 10 Gyr) along with the old-HB stars and RR 
Lyrae variables. 

5. Summary and conclusions 

We have presented a large area study of the field population 
in the Fornax dwarf spheroidal galaxy, based on BVI data for 
about 40000 stars. The size of our sample, together with the 
wide photometric baseline employed in this work, provide new 
information on the stellar content of Fornax. 

One of the most distinctive features in our diagrams is 
the conspicuous young main sequence. In this paper we have 
shown that the plume of stars just above the red clump is made 
up of intermediate mass stars (2.4-2.9 Mq) burning helium 
in the core, counterparts of the young main sequence stars in 
the age range 0.3-0.4 Gyr. The comparison with isochrones 
suggests us that these blue-loop stars must be as metal-rich as 
[Fe/H]^ —0.7, which represents an important constraint for 
the metal enrichment history in Fornax. 

An extended upper AGB tail and a prominent red HB clump 
testify the presence of a dominant intermediate-age population 
in the age range 2-10 Gyr, corresponding to 0.9-1.4 M Q stars. 
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From the difference in the mean (V — I) colors of the red clump 
and the RGB at the same luminosity, we have estimated a mean 
age 5.4 ±1.7 for the bulk of the intermediate-age population in 
Fornax, in agreement with the morphology of the MS turnoffs 
in WFPC2 color-magnitude diagrams (Buonanno et al. 1999). 



This suggests that the location of the red HB clump may indeed 
prove to be a useful age indicator for distant LG galaxies. 

About 0.2 mag below the red clump, an extended HB is 
indicative of an old population. In particular, our data point to 
the presence of blue HB stars in the outer regions. Together 
with previous detection of RR Lyrae, this provides evidence 
for a minority field population that is as old and metal-poor as 
that in the Fornax globular clusters. The Fornax dSph clearly 
started forming stars in a halo nearly at the same epoch when 
most of its surrounding clusters were formed. 

Evolutionary phases that gave barely discernible features 
in small field observations are easily measurable in our color- 
magnitude diagrams. We could reliably measure the AGB 
bump, a small clump produced by a clustering of stars at the 
base of the AGB, at My — —0.4. Measurements of such minor 
evolutionary features may provide useful tests of stellar models 
for stars of different masses and metallicities. 

The sharp cutoff in the luminosity function of Fornax has 
been used to estimate its distance using the RGB tip method. 
The corrected distance modulus of Fornax, (m — M)o = 
20.70 ± 0.12, agrees with previous determinations. This esti- 
mate is confirmed by the mean level of old horizontal-branch 
stars. By measuring the average magnitude of the red HB (dis- 
tinct from the red clump) we estimated a distance modulus 
(m — M)q = 20.76 ± 0.12 on the distance scale of Lee et 



al. ( |1990| ). 

Fornax, as many other dSph, has been known for a long 
time to have a wide RGB color distribution. The "color scat- 
ter" has been usually taken to represent an abundance spread. 
We have analyzed in detail the color distribution of the red giant 
stars across the fiducial line, and found that it is reasonably well 
fitted by a two-component model. This approximately bimodal 
distribution is remarkably similar in all fields. About 70% of 
the red giants belong to an intermediate-age RGB component 
which is itself wider than expected from instrumental errors. 
By comparing the bulk of the Fornax RGB with the ridge lines 
of standard globular clusters, we have estimated a mean metal- 
licity [Fe/H]= —1.39 ± 0.15. This nominal value should be 
corrected for the age difference between the Fornax popula- 
tion and the Milky Way globulars. Accounting for an age dif- 
ference of 10 Gyr, we find an age-corrected mean metallic - 
ity [Fe/H]= —1.0 ± 0.15 for the dominant intermediate-age 
population of Fornax. Interestingly, this is also the metallic - 
ity found for Sagittarius, the nearest Milky Way dSph satellite 
that has luminosity comparable to that of Fornax. The intrinsic 
color scatter of stars in the RGB main component is <J(b-i) = 
0.06±0.01 mag implying a relatively modest metallicity spread 
crpe/H] = 0.12 ± 0.02 dex. Then the (2a) metallicity range for 
the bulk population of Fornax is -1.25 < [Fe/H] < -0.75 if a 
correction for the mean age is applied. The secondary compo- 



nent or "bluer tail" is also quite broad. In principle, these bluer 
stars could be either young or old and metal-poor. 

Star counts of different subpopulations at various locations 
confirm and extend the evidence for radial population gradients 
emerged in previous studies. Recent star formation is clearly 
concentrated in the central regions, though with some degree 
of asymmetry (e.g., SHS98). Old stars are more easily seen in 
the outer fields. A blue HB population can be noticed in our 
outermost field, coming from the minority old, metal-poor field 
component. The stars populating the blue side of the wide RGB 
closely follow the spatial distribution of the old-HB stars. This 
is perhaps our most important finding, since it demonstrates 
that the bluer RGB stars are themselves old and metal- poor, 
and clearly establishes the nature of the wide RGB of Fornax. 
Thus the roughly bimodal color distribution can be interpreted 
as a metallicity distribution, implying that the bulk of the For- 
nax galaxy was built during two rather distinct star-forming 
epochs. The older population has [Fe/H]?s —1.8 dex with 
(±2er) and an wide abundance range —2.2 < [Fe/H] < —1.4. 

The emerging picture is one in which the evolution of For- 
nax is characterized by two major star formation epochs, each 
consisting of many episodes. The first episode took place at an 
early epoch, being presumably coeval to the birth of the old 
galactic globular clusters, from metal-poor gas. After a rela- 
tively quiescent period, Fornax formed the bulk of stellar pop- 
ulations between 7 and 2.5 Gyr ago from the pre-enriched gas. 
Star formation continued at a lower rate in the central regions 
until as recently as 10 8 yr ago. The modest internal abundance 
spread found in each main population seen in the metallicity 
distribution, and the different metallicities of populations of 
different age, trace the progressive metal enrichment and rep- 
resent the basis for an age-metallicity relation in Fornax. The 
constraints found in this paper provide the physical input for a 
quantitative analysis of the star formation and chemical enrich- 
ment history of Fornax, which will be done in a forthcoming 
study using the methods of stellar population synthesis. 
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